Switching and Jamming Transistor Effect for Vortex Matter in Honeycomb Pinning 

Arrays with ac Drives 
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We show that a remarkable variety of dynamical phenomena, including switching, polarization, 
symmetry locking, and dynamically induced pinning, can occur for vortices in type-II superconduc- 
tors in the presence of a honeycomb pinning array and an ac or combined ac and dc drive. These 
effects occur at the second matching field where there are two vortices per pinning site, and arise 
due to the formation of vortex dimer states in the interstitial regions of the honeycomb array. The 
orientation of the pinned and moving vortex dimers can be controlled externally by the application 
of a drive. We term this a polarization effect and demonstrate that it can lock or unlock the vortex 
motion into different symmetry directions of the underlying pinning lattice. If the moving vortices 
are locked into one direction, the motion can be switched into a different direction by applying an 
additional bias drive, producing sharp jumps in the transverse and longitudinal velocities. Further, 
the dc vortex motion in one direction can be controlled directly by application of a force in the 
perpendicular direction. When the moving dimers reorient, we find a remarkable dynamical pinning 
effect in which the dimers jam when they become perpendicular to the easy flow direction of the 
pinning lattice. Since application of an external field can be used to switch off the vortex flow, we 
term this a jamming transistor effect. These effects do not occur in triangular pinning arrays due 
to the lack of the n-merization of the vortices in this case. The switching and dynamical pinning 
effects demonstrated here may be useful for the creation of new types of fluxtronic devices. 

PACS numbers: 74.25.Qt 



I. INTRODUCTION 



The dynamics of vortices in type-II superconductors 
has attracted enormous interest since they can provide in- 
sight into nonequilibrium phases of driven particles mov- 
ing over a quenched landscape and properties of vortex 
dynamics that can be useful for creating new types of 
fluxtronic devices. In terms of basic science, driven vor- 
tices are an ideal system in which to study different types 
of nonequilibrium phase transitions P, 0, [3, H, H, 0, Q ■ It 
appears that nonequilibrium systems can undergo tran- 
sitions between specific types of driven phases similar to 
the transitions among equilibrium phases that occur as 
a parameter is tuned. Vortices driven over a randomly 
disordered substrate can pass through a pinned phase, a 
disordered or plastic flow phase with large fluctuations 
[Hi S Q J and different types of moving phases which may 
be partially ordered in one direction to form a moving 
smectic [l,[a,[3l or have moving crystal symmetry [3, 4, 7]. 
Other types of phases are also possible within the plastic 
flow regime, including chaotic Q, reversible, and irre- 
versible flows 0. For driven vortices in samples with 
periodic pinning arrays, an exceptionally rich variety of 
distinct dynamical behaviors a ppear [TqI . [Ill , [l^ . [T3l . [H, 
m, M [13, M M M M M MM, llI~Periodic pinning 
arrays in superconductors can be created with nanoholes 
dE M, MMMM, M M [iy, iSS,, M, M M or mag- 
netic dots [3^, Eo, EH . Typically, strong matching effects 
are observed in these samples as anomalies in the crit- 
ical current or other transport measurements when the 
number of vortices equals an integer multiple or rational 
fraction of the number of pinning sites. At the matching 



field B^, there are an equal number of vortices and pins. 
Above the first matching field, the additional vortices cor- 
responding to the excess density B — B^ can be located 
in the interstitial regions between the occupied pinning 
sites already occupied by vortices [13, HI, [30l. l33l. [34 | 
or multiple vortices can be trapped at each pinning site 
(29I [30I m, HI]. If interstitial vortices are present, 
they are more mobile than the pinned vortices under an 
applied drive and have a lower depinning threshold, so 
the interstitial vortices can move around the pinned vor- 
tices as observed in experiments [ij, [l^ and simulations 
[Tol . [l3j . The interstitial vortex motion can be probed 
by adding an ac drive to a dc drive to produce Shapiro 
steps in the current- voltage curves 15| . The symmetry of 



the underlying pinning array can lock the vortex motion 
into particular directions when the vortices are driven 
at various angles with respect to the pinning lattice, as 
predicted in simulations [l^ and demonstrated in exper- 
iment [13, Hll, In addition to the superconducting 
system, vortices interacting with periodic substrates can 
also be realized in Bose-Einstein condensates with a peri- 
odic optical trap array [i^, [H, , opening a new avenue 
for the study of vortex matter on dynamical substrates. 

Studies of vortices interacting with periodic pinning ar- 
rays have primarily concentrated on triangular or square 
arrays; however, other types of pinning geometries such 
as honeycomb or kag ome arrays can also be realized 
in experiment [38l . |41| |. A honeycomb pinning array is 
constructed from a triangular array by removing 1/3 of 
the pinning sites, while a kagome array is generated by 
removing 1/4 of the pinning sites. Experiments have 
shown that honeycomb and kagome pinning arrays ex- 
hibit anomalies in the critical currents at non-matching 
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fields as well as at matching fields [H, |4l|. In honey- 
comb pinning arrays, non-matching anomalies occur at 
fields B/Bcji = n + 1/2, where n is an integer. This sug- 
gests that a portion of the vortices are trapped in the 
large interstitial regions in such a way as to match the 
original undiluted triangular pinning array. Simulations 
also reveal matching effects at integer and half integer 
matching fields up to B/B^ — 5 for honeycomb arrays, 
and indicate that this type of matching can occur when 
multiple interstitial vortices are trapped at the large in- 
terstitial sites in the honeycomb array and form either 
ordered vortex molecular crystals or vortex plastic crys- 
tals (4^]. In vortex molecular crystals, the vortices in the 
large interstitial site behave as n-mer objects that have 
a director field which can be oriented in different direc- 
tions with respect to the pinning lattice. The n-mers can 
have long range order and may experience ferromagnetic 
coupling that aligns the n-mers in a single direction, an- 
tiferromagnetic coupling that causes neighboring n-mers 
to orient perpendicularly to each other, or other types of 
coupling that produce ordered phases [43| • Vortex molec- 
ular crystals form at integer and half integer matching 
fields when B / B^ > 2 for honeycomb pinning arrays and 
at inte ger multiples of B/B^ = 1/3 for kagome pinning 
arrays [43]. Vortex plastic crystals occur when the tem- 
perature is high enough to destroy the orientational order 
of the vortex n-mers. Experiments show that as the tem- 
perature increases, certain commensuration peaks disap- 
pear, which is consistent with the formation of plastic 
vortex crystals [38, 45]. 

The vortex n-merization in honeycomb pinning arrays 
produces a novel dynamical symmetry breaking effect 
[461 lirf . At B/B^ = 2, the vortices form dimers with 
a broken ground state symmetry. When a longitudinal 
drive is applied, there is a positive or negative trans- 
verse velocity response from the moving dimers depend- 
ing on the direction of the original symmetry breaking. 
At incommensurate fields near B/B^ = 2 where dimers 
are present, the ground state is not globally symmetry 
broken; however, under a longitudinal drive the moving 
dimers can dynamically organize into one of two possible 
flow states with a transverse response. When the longitu- 
dinal drive is held fixed at a constant value, the direction 
of the transverse flow can be switched with an externally 
applied dc transverse drive. Both the transverse and lon- 
gitudinal velocities change abruptly when the flow direc- 
tion switches. If the amplitude of the transverse drive 
is too weak, the system remains locked in the original 
transverse flow direction. When an ac transverse drive is 
applied in addition to a flxed dc longitudinal drive, the 
transverse portion of the vortex flow oscillates, resulting 
in an ampliflcation of the transverse response in certain 
regimes. In these regimes, an increase of the transverse 
ac force can cause the longitudinal vortex velocity to de- 
crease. The pronounced switching and amplification be- 
haviors of the vortex dimers may be useful for producing 
transistor like devices [H, |4§] . 

In this work, we study the vortex dimer response in a 



sample with a honeycomb pinning lattice subjected to 
ac or combined ac and dc drives. We flnd switching 
and jamming phenomena and a rich variety of dynam- 
ical behaviors. Many of our results should be relevant 
for systems other than vortices on periodic pinning ar- 
rays. Simulations, theory, and experiments have shown 
that colloidal particles interacting with a periodic sub- 
strate can exhibit an n-merization effect at fillings where 
there is more than one colloid per substrate minimum 
[sol [sil m, [si] . These states have been termed colloidal 
molecular crystals and have man y of the same properties 
as the vortex molecular crystals [isl l46l . |47| . Our results 
may also apply to other systems such as molecular dimers 
driven over periodic substrates, friction systems [13], or 
charged balls on periodic substrates [s^. The symmetry 
locking effect predicted for vortices moving over periodic 
arrays flSJ] has been observed in experiments with col- 
loids moving over periodic arrays [56|. We thus expect 
that the results described in this work should be generic 
for colloids on honeycomb arrays. 



II. SIMULATION 

We simulate a two-dimensional system of size L x L 
with periodic boundary conditions in the x and y direc- 
tions containing 7V„ vortices at a density n„ = N^/Lp'. 
The vortex density n^ is proportional to the magnetic 
field B = Bi. Vortices can move under the influence of 
the Lorentz force from an applied current J or due to 
interactions with the other vortices. In a type-II super- 
conductor the vortex core is not superconducting, so the 
vortex motion is overdamped due to dissipation of en- 
ergy by the core. The motion of vortex i at position 
is governed by the following equation of motion: 

^F^+Ff + F'^^ + F''^ + Ff. (1) 

Here the damping constant rj — (f)Qd/2Tr£^^pj\j, where d 
is the sample thickness, (f)o — h/2e is the flux quantum, 
^ is the superconducting coherence length, and pat is the 
normal state resistivity of the material. The repulsive 
vortex-vortex interaction force is 

Fr = £/oi^i(^)R.„ (2) 

where Ki is the modified Bessel function, A is the Lon- 
don penetration depth, /o ~ (/)Q/(27r/ioA'^), the distance 
between vortices i and j is Rij — ]Ri — Rj ], and R^j — 
(Ri — Rj)/i?.y. We measure lengths in units of A and 
take L — 24A so that the system size is 24A x 24A. Since 
the vortex-vortex interaction force falls off rapidly with 
Rij, we impose a long range cutoff at 7A. For very short 
Rij a cutoff is placed at 0.1 A to avoid a divergence in the 
force. 
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FIG. 1: Vortex positions (filled dots) and pinning site lo- 
cations (open circles) for a honeycomb pinning array with 
B/B^ — 2.0. Here the interstitial vortices form a dimer 
state. A circular ac drive is applied with A = 0.075 and 
T = 3.14X 10^ simulation time steps, (a) In the initial state at 
t — the dimers are aligned along the y-direction ed, 6 — -k 12. 
(b) At tjT = 0.25, the dimers are tilted to 6* = 57r/6. (c) At 
tjT = 0.54, the dimers are in the process of switching to a 
new symmetry direction, (d) At tjT = 0.6, the dimers have 
switched to the orientation 6 = ln/Q. 



The force from the vortex-pinning site interactions is 
F^. We place Np pinning sites, with a pinning density 
of Up = Npjl? ^ in a honeycomb array formed by remov- 
ing every third pinning site from a triangular array. The 
pinning sites have a maximum force of Fp — 0.85 and a 
fixed radius of — 0.35A, which is small enough that 
only one vortex is trapped per pinning site. Each pin- 
ning potential is modeled as an attractive parabolic well, 
giving a force 



'P T? nip) f 



R 



'ik 



Here B is the Heaviside step function, Rl is the loca- 



(3) 



(p) 

ik 



IR^-Ri^^l, and R^^^ 



tion of pinning site fc, R, 

(R, - R^f^-Rifc^- The matching field is defined as 
the field with iV„ = Np where there is one vortex per 
pinning site. We fix = 0.3125(^o/A^. In Fig. 1 we 
illustrate the pinning site positions and vortex locations 
for a honeycomb pinning array at B / B^ = 2.0. 

The thermal force F^ is modeled as random Langevin 
kicks with the correlations (Ff ) = and (Ff (t)Fj(i')) = 
2r]kBT5{t — t')5ij. We obtain the initial vortex configura- 



tions by simulated annealing. In this study we anneal by 
starting at Jt = 3.0 and spend 5000 simulation time steps 
at each temperature before reducing the temperature in 
increments of 5 fx — 0.002. Slower annealing times do 
not change the final T — Q vortex configurations. 

Once the vortex configuration is initialized, we apply 
a uniform dc drive F'^'^ = F'^'^y and ac drive F"'^. Unless 
otherwise noted, the ac force has the form 



F°= = Afo [sin(27rt/T)A + cos(27rt/r)y] 



(4) 



with amplitude A and period T. We measure the average 



vortex velocity response Vx = N^^ (I2f=i ' ^) ^^'^ 



^v^(I2f=i'^i • y), where is the velocity of vortex i. 
We note that the ability to apply currents in both x and 
y directions and measure the voltage response has been 
demonstrated experimentally |2l|, . 



III. EXTERNAL CONTROL OF THE VORTEX 
DIMER ORIENTATION 

We first consider a honeycomb pinning lattice at 
B/ Bcf, — 2.0 with a circular ac drive and no dc drive. At 
this matching field, the interstitial vortices form dimer 
states in the large interstitial regions of the honeycomb 
pinning lattice, as shown in Fig. 1(a) and studied previ- 
ously in Ref. [45[ . The dimers have an effective ferromag- 
netic coupling and are all aligned in the same direction. 
In previous work we showed that in the absence of a 
drive, the ground state configuration of the dimers can 
be aligned in one of three directions Here we show 
that an external ac drive can couple to the dimers and 
change their orientation provided that the ac amplitude 
is small enough that the interstitial vortices do not depin. 

In Fig. 1 we illustrate the time evolution of a system 
under an ac drive with A = 0.075 and T = 3.14x 10® sim- 
ulation time steps. The dimers are initially aligned in the 
?/-direction as shown in Fig. 1(a). We denote the direction 
of the dimer orientation as 9, measured on the unit circle. 
Here, the initial orientation \s 9 = 7r/2. The ac drives in- 
duces two types of dimer motion. The first is a rotation 
of the dimers in phase with the ac drive and the second 
is an abrupt switching where the dimers suddenly change 
their orientation to a new symmetry direction. We char- 
acterize these different types of motion by examining the 
vortex positions and trajectories as well as the time traces 
of Vx and Vy. Fig. 1(b) shows a snapshot of the dimer po- 
sitions at t/T = 0.25 when the dimers have tilted to align 
aX 9 = 57r/6. At t/T = 0.54 the dimers are in the process 
of switching between 6* = 57r/6 and 9 = 7tt/6, and align 
as shown in Fig. 1(c) to form an overall stripelike pat- 
tern that incorporates the pinned vortices. In Fig. 1(d) 
at t/T = 0.6, the dimers have switched to an alignment 
oi 9 = 7tt/6. The time sequence shows that the dimers 
rotate in phase with the circular ac drive. 

In Fig. 2 we plot the vortex trajectories corresponding 
to the system shown in Fig. 1. The first two switches 
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FIG. 2: Vortex positions (filled dots), pinning site locations 
(open circles), and vortex trajectories (lines) for the same 
system in Fig. 1. (a) The trajectories during a time interval 
of 0.6r when two switching events occur, (b) The trajectories 
during a time interval of I.IT showing the complete vortex 
orbits. The interstitial vortices follow the same orbits for all 
later times. 



in the cycle (during a time interval of 0.6T) appear in 
Fig. 2(a), while Fig. 2(b) indicates that the dimers move 
in a closed triangular orbit during a complete cycle. A 
combination of smooth and switching dimer motion ap- 
pears in Vx and Vy as a function of time in Fig. 3. The 
sinusoidal envelopes of the velocity traces are interlaced 
with spike features formed when the dimers abruptly 
switch to a new orientation. The first jump at t/T = 0.2 
occurs at the transition from the 6 = tt/2 orientation 
shown in Fig. 1(a) to the 9 — 5tt/6 orientation shown 
in Fig. 1(b). The lower vortex in each dimer shifts in 
the positive x direction with a smaller displacement in 
the negative y direction, while the other vortex in the 
dimer shifts in the negative y direction with a smaller dis- 
placement in the negative x direction. The sum of these 
motions gives a positive Vx pulse and a simultaneous neg- 
ative Vy pulse in Fig. 3. In the next switching event at 
t/T — 0.54, the upper vortex in each dimer shifts in the 
negative y direction while the other dimer vortex shifts 
in the positive y direction, reducing the net value of Vy. 
At the same time, both vortices in the dimer move in the 
negative x direction, enhancing the magnitude of Vx as 
shown in Fig. 3. This process repeats each period. The 
switching behavior occurs since the vortex dimers sit in 
a caging potential with a sixfold modulation created by 
the repulsive interactions with the six vortices located at 
the neighboring pinning sites. This modulation provides 
an energy barrier for the switching of the dimer from one 
orientation to another. The field-induced switching effect 
we observe is similar to liquid crystal behavior where the 
orientation of rod shaped molecules can be altered with 
external fields. Thus, we term the dimer switching a po- 
larization effect, and note that it could be used to create 
numerous novel vortex devices. For example, the orien- 
tation of the dimers can be rotated in such way as to 
maximize the critical current in certain directions. Addi- 
tionally, logic states could be assigned with different logic 
values corresponding to different dimer orientations. 




FIG. 3: The instantaneous vortex velocity vs time for the 
system in Fig. 1. (a) Vx- (b) Vy. The spikes correspond to the 
abrupt switching events where the dimers change orientation. 



The pinned switching behavior shown in Fig. 3 only 
occur under certain conditions. If the ac amplitude A is 
small, the dimers are unable to switch and remain locked 
in a single orientation. Conversely, if the ac amplitude is 
too large, the dimers escape from the interstitial caging 
potential and channel along one of six angles: 6 — 7r/6, 
tt/2, 57r/6, 7n/6, 37r/2, or ll7r/6. In Fig. 4(a,b) we plot 
Vx and Vy versus time for the same system in Fig. 3 
but with a smaller ac amplitude of A = 0.03 which is too 
weak to induce switching. The velocity curves are smooth 
and merely follow the ac input drive as the vortices shift 
slightly in the interstitial sites. When the external drive 
is large enough, the dimers depin from the large intersti- 
tial regions and channel along one of the six easy-flow di- 
rections of the lattice, giving the velocity response shown 
in Fig. 4(c,d) for A — 0.125. The channeling effect is very 
pronounced due to the significant barriers that the vor- 
tices must overcome in order to lock their motion into a 
new easy-flow direction. In fact, the locking is so strong 
that when the vortices switch to a different flow direction 
they skip over every other easy-flow direction, so that 
only three of the possible six symmetry flow directions 
are followed in a given realization. Which of the two sets 
of three flow directions that the system chooses depends 
on the initial conditions. In Fig. 4(c,d) at ^ = 0.125, the 
dimers initially flow along 6 — 7r/2, similar to the motion 
shown in Fig. 5(a). Figure 4(c,d) indicates that during 
the first portion of the drive cycle, —0.15 < t/T < 0.15, 
Vy is finite and Vx is zero, as expected for motion that 
is locked in the y direction. For 0.15 < t/T < 0.5, the 
X component of the ac drive becomes strong enough to 
reorient the dimers so that they channel along a differ- 
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FIG. 4: (a) Vx and (b) Vy versus time for the same system in 
Fig. 1 with A = 0.03 where the dimers do not switch, (c) Vx 
and (d) Vy for the same system with A = 0.125 where the ac 
amplitude is high enough to induce a sliding of the dimers. In 
each driving period, the dimers initially slide along 6 = ir /2 in 
the positive j/-direction, become pinned, and then slide along 
6) = ll7r/6 and e = 7tv/6 in turn. 



ent symmetry direction. The vortices lock to 6 = ll7r/6 
and have the same type of motion illustrated in Fig. 5(e). 
During 0.5 < t/T < 0.85, the vortex motion switches and 
flows along 6 = 7tt/6. Finally the vortex motion returns 
to 6* = 7r/2. In this cycle, the vortices never flow along 
6 = n/G, 37r/2, or 57r/6. By changing the initial condi- 
tions slightly, we can obtain vortex motion along those 
three angles only (and not along tt/2, ll7r/6, or 7tt/6) 
with a fifty percent probability. For A — 0.125, there 
are small windows of pinned phases between the sliding 
locked phases. In Fig. 6 we plot the vortex trajecto- 
ries during an entire drive period showing the triangular 
pattern formed as the vortices flow along three different 
directions. 

For ac amplitudes 0.11 < A < 0.2, Fig. 7(a) shows that 
only three of the six possible flow directions appear and 
that Vx and Vy simultaneously drop to zero between the 
sliding regimes in each period, indicating the existence 
of a temporary pinned phase. For A > 0.2 the system 
passes directly from one siding state to another with- 
out pinning, as indicated in Fig. 7(b), and flow occurs 
in all six possible directions. Figure 7(b) shows that at 
A = 0.3, Vx and Vy have discontinuities when the system 
switches from one flow direction to another, but that Vx 
and Vy are never simultaneously zero. Within a given 
flowing locked phase, Vx and Vy are not fixed but can 
show considerable variations. 

For 0.2 < A < 0.335, only the interstitial vortices are 
moving, while for A > 0.335 a portion of the vortices lo- 
cated at the pinning sites become mobile during certain 
parts of the ac drive cycle. This produces additional dis- 
ordered flow phases between the flowing locked phases. 
We plot Vx and Vy during one drive cycle at ^ = 0.35 
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FIG. 5: Vortex positions (filled dots), pinning site locations 
(open circles), and vortex trajectories (lines) for the system 
in Fig. 1 at A = 0.35 during half of a drive cycle. The points 
from which the images are taken are marked with bold letters 
in Fig. 8(a). (a) The dimers initially move in the j/-direction, 
9 — n/2. (b) Transition to the disordered flow state in which 
a portion of the pinned vortices depin. (c) The dimers lock 
to 5 = 7r/6. (d) Transition to another disordered flow state 
where a portion of the pinned vortices depin; the average 
velocity is in the x direction, (e) The dimers lock to 6 = 
ll7r/6. (f) Transition to a disordered flow state in which a 
portion of the pinned vortices depin. The remainder of the 
drive cycle has the same appearance as (a-f) but with the 
motion of the vortices reversed to follow — 3tt/2, 7n/6, and 
57r/6. 



in Fig. 8(a) and show the corresponding vortex trajecto- 
ries in Fig. 5. At the beginning of the drive cycle, the 
dimers are locked to flow along 9 = 7r/2, as shown in 
Fig. 5(a). The system does not jump directly into an or- 
dered channel flow of dimers along 6 — tt/6, but instead 
first passes through an intermediate flow state near point 
b in Fig. 8(a). The vortex trajectories are disordered, 
as shown in Fig. 5(b). Within this regime a portion of 
the vortices at the pinning sites depin due to the com- 
bination of the vortex-vortex interaction forces and the 
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FIG. 6: Vortex positions (filled dots), pinning site locations 
(open circles), and vortex trajectories (lines) for the system in 
Fig. 4(c,d) at A = 0.125 during a complete drive cycle. The 
dimers slide along three different angles so the trajectories 
form a triangular pattern. 

applied drive. The result is a random flow which ap- 
pears as enhanced fluctuations in Vx and Vy. From the 
random phase, the system reorders into a locked flow- 
ing state along 6 = tt/6 illustrated in Fig. 5(c), where 
all the pinning sites are again occupied and the dimers 
move in an ordered fashion. The system then enters an- 
other disordered flow regime shown in Fig. 5(d). The 
transition into this disordered phase produces a jump in 
Vx- As the ac drive period progresses, the vortices lock 
into another flowing state along 6 = IItt/G illustrated in 
Fig. 5(e), followed by the disordered flow phase shown 
in Fig. 5(f). During the second half of the drive cycle, a 
similar series of ordered and disordered phases occur as 
the vortices successively lock into 9 = 37r/2, 9 — 7tt/6, 
and 9 = 57r/6, before finally returning back to the start- 
ing point of locked flow along 9 ~ tt/2. In all cases, the 
random flow phases are associated with increased fluctu- 
ations of Vx and Vy. The vortex velocity distributions 
are broadened in the random flow phase, whereas in the 
ordered flow phases the velocities are synchronized. 

In Fig. 8(b) we plot Vx and Vy versus time for A = 
0.405 where only a small window of each locked slid- 
ing phase occurs between the random flow regimes. For 
A > 0.42, the locked sliding phases vanish and the flow is 
disordered throughout the ac drive cycle, although some 
slight channeling persists within the random flow phase. 

We summarize the results for < A < 0.45 in a dy- 
namic phase diagram in Fig. 9, where we plot the phases 
for A versus time during half an ac drive period. For 




FIG. 7: Vx (heavy line) and Vy (light line) versus time for 
the system in Fig. 1. (a) At A = 0.16 there are three pinned 
phases between the three sliding phases in each period, (b) 
At A = 0.3, the dimers are always sliding along the different 
symmetry directions and there are no pinned phases. 



A < 0.05 we find a pinned locked (PL) phase in which all 
the vortices remain pinned and the dimers do not show 
any switching behavior. The dashed line at A = 0.05 
separates the PL from the pinned switching (PS) regime 
where the dimers remain pinned but change orientation. 
For 0.11 < A < 0.16 the PS phase is lost over much 
of the drive cycle and is replaced by a series of sliding 
states along every other symmetry angle. Here, the vor- 
tices follow 9 = Tr/2 and 9 = llir/G in this half of the 
drive cycle, separated by the remnants of the PS phase. 
The PS regime disappears completely for A > 0.16 and 
for 0.16 < A < 0.3, the vortices lock into each symme- 
try angle in turn, passing from = 7r/2to0 = 7r/6to 
9 = ll7r/6, and so forth. For 0.325 < A < 0.435, the 
sliding phases are separated by intrusions of the random 
flow (R) phase which increase in extent until the sliding 
states are completely lost for A > 0.435. 



A. High ac Amplitude and Order-Disorder 
Transitions 

The vortices depin and flow in a random phase when 
the ac amplitude is large enough to overcome the pinning 
force. For even higher ac drive amplitudes, we flnd that 
the system enters a new dynamical regime. The vortex 
flow changes from the channeling of interstitial vortices 
between occupied pinning sites to a regime where all the 
vortices take part in the motion and a portion of the vor- 
tices move in locked channels that pass through the pin- 
ning sites. For ac amplitudes just above A = 0.435 where 
the dimer channeling is lost, the system is disordered and 
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FIG. 8: Vx (heavy lines) and Vy (light lines) versus time for 
the same system in Fig. 1. (a) At A — 0.35, a portion of 
the vortices at the pinning sites start to depin and generate 
random phases between the locked sliding phases. The labels 

a, b, c, d, e, and f correspond to the portions of the period 
illustrated in Fig. 5. The disordered flow phases appear at 

b, d, and f during the transitions between different locked 
sliding states, (b) At A = 0.405, the random flow phases are 
broader. 
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FIG. 9: The dynamic phase diagram for the system in Fig. 1 
showing ac amplitude A vs time. For low A the dimers are 
in the pinned locked (PL) state. For 0.05 < A < 0.11 the 
dimers are in the pinned switching (PS) state where they 
remain pinned but switch between different orientations as 
shown in Fig. 2 and Fig. 3. The dashed line separates the 
PL and PS phases. For 0.11 < A < 0.16, the vortices pass 
through a series of sliding phases along 9 = it/2, 6 — it/6, 
and 6 = ll7r/6. The sliding phases are separated by pinned 
phases for 0.11 < A < 0.16 and by random flow (R) phases 
for 0.325 <A< 0.435. 



the vortices undergo plastic flow in which some vortices 
are temporarily trapped at the pinning sites. The fea- 
tures in the velocity versus time curves are washed out 
and we find only a weak channeling of the vortex motion. 
In this regime, the vortex lattice is heavily defected. As 
the ac amplitude is further increased, the effectiveness of 
the pinning diminishes and the vortices begin to partially 
reorder along certain symmetry directions of the pinning 
lattice. 

To quantify the reordering we use a Voronoi construc- 
tion to measure the fraction of six-fold coordinated vor- 
tices Pa — X]il="i ^i'^i ~ 6)' where Zi is the coordina- 
tion number of vortex i. For a perfectly triangular vortex 
lattice, Pq ~ 1. In Fig. 10(a) we plot Vx and Vy vs time 
for A = 0.43, where the vortices are in the random flow 
(R) regime. Small features in Vx and Vy indicate the 
persistence of partial channeling even in the random flow 
phase. Figure 10(b) shows the corresponding measure- 
ment of Pq, which has an average value of Pq « 0.5 due 
to the highly defective nature of the vortex lattice in the 
random flow phase. At ^ = 0.5, Fig. 10(c) indicates that 
the channeling features disappear and are replaced by 
new steplike features in Vx and Vy. Figure 10(d) shows 
that the steps are correlated with six maxima in Pg, in- 
dicating that the vortices are partially ordered along the 
steps and that a different type of channeling effect is oc- 
curring. The overall vortex lattice remains disordered 
throughout the drive cycle; however, when the drive falls 
on a step, Pq ~ 0.68, while between the steps, Pq ~ 0.5. 

At ^ = 0.75, shown in Fig. 11(a), the steps in Vx and 
Vy are much denser and are associated with a series of 
peaks and dips in Pq, shown in Fig. 11(b). Between the 
steps, Pg « 1, indicating a triangular vortex lattice, while 
on the steps, the vortices are less ordered and Pg « 0.68. 
The changes in Pg indicate that the vortices undergo a 
series of order-disorder transitions during the ac drive 
cycle. The average minimum amount of vortex order 
present at A = 0.75 is larger than average maximum 
vortex order ai A — 0.5, as indicated by comparing Pg in 
Fig. 11(b) and Fig. 10(d). In Fig. 12(a) we show a blowup 
of Vx and Vy for 0.05 < t/T < 0.2 from Fig. 11(a) to high- 
light the order-disorder transition that occurs when the 
vortices enter and exit a single velocity step. We plot the 
corresponding Pg values in Fig. 12(b). Within the non- 
step intervals, where Vx and Vy are smoothly increasing 
or decreasing, Pg « 1, indicating that the vortices are 
mostly sixfold coordinated. The dashed lines in Fig. 12 
outline the step located at 0.072 < t/T < 0.102. A dip in 
Pg occurs as the system both enters and exits this step, 
indicating enhanced disorder in the vortex lattice. On 
the velocity step, Pg increases to Pg w 0.85 but the vor- 
tex lattice structure is not completely ordered. Outside 
of the step, the vortices regain a nearly triangular order- 
ing until reaching the next velocity step at t/T « 0.15. In 
Fig. 13 we plot the vortex trajectories for the system in 
Fig. 12 before, during, and after the step, where the step 
onset tstep/T = 0.072. Below the step. Fig. 13(a) shows 
that the vortex trajectories are disordered, while on the 
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FIG. 10: (a) Vx (dark line) and Vy (light line) versus time 
for the system in Fig. 9 sd, A = 0.43, where the flow is dis- 
ordered. Although the system is in the random phase, there 
is still some mild channeling of the vortices which appears as 
modulations of Vx and Vy. (b) The corresponding fraction 
of six-fold coordinated vortices Pe vs time. For a completely 
triangular lattice Pg = 1.0. Here Pe ~ 0.5, indicating that 
the vortex lattice is disordered throughout the whole driving 
cycle, (c) Vx (dark line) and Vy (light line) versus time for 
the system in Fig. 9 ai A — 0.5. The channeling effects are 
completely destroyed; however, new steplike features are be- 
ginning to appear, (d) The corresponding Pe vs time shows 
enhancements that are correlated with the steplike features 
in Vx and Vy, indicating that the vortices are gaining some 
partial ordering on the steps. 
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t/T 

FIG. 11: (a) Vx (dark line) and Vy (light line) vs time for 
the same system in Fig. 10 with a higher ac amplitude of 
A — 0.75. (b) The corresponding Pe vs time shows a series of 
order to disorder transitions. 



step Fig. 13(b) shows more ordered vortex motion with 
a partial locking to 9 = 7r/3. The vortices channel along 
the pinning sites rather than between them. Above the 
step, the vortex motion is again disordered, as illustrated 
in Fig. 13(c). 

The other velocity steps in Fig. 11(a) also correspond 
to regimes where the vortices channel along the pinning 
sites at different angles. This type of channeling motion 
along pinning sites was studied previously for vortices in 
square pinning arrays at fields close to the matching field 
B/Bcf, = 1.0 In that system, a fixed dc drive was 
applied in the longitudinal direction and a transverse dc 
drive Fj was gradually increased. A series of steps in 
the transverse velocity appeared, forming a devil's stair- 
case structure. The velocity steps form when the vortex 
motion is locked to certain angles with respect to the 
symmetry axis of the pinning array. Specifically, for a 
square array steps form at drive ratios of F^/F^ = m/n, 
where m and n are integers. For a triangular pinning ar- 
ray, velocity steps occur when Ff / F^ = ^Jrri/2n -f 1. For 
the honeycomb array, we expect a similar set of velocity 
steps to occur at the same drive ratio values as in the tri- 
angular array, but unlike the triangular array, the steps 
should be particularly pronounced when the drive angle 
is along = 7r/3, = 7r/2, 6 = 27r/3, or the complements 
of these angles. In Ref. [l^, a series of order-disorder 
transitions were observed that are similar to the order- 
disorder transitions we find in the honeycomb array. 

In Fig. 14(a) we plot the Voronoi construction of a 
snapshot of the vortex configuration at (tstep — t)/T ~ 
0.06, below the step highlighted in Fig. 12. In the non- 
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FIG. 12; A blowup of the region 0.05 < t/T < 0.2 from Fig. 11 
with A — 0.75 highlighting the order to disorder transitions 
across a step, (a) Vx (dark line) and Vy (light line) versus 
time, (b) The corresponding Pg vs time. The dashed lines 
highlight the transitions into and out of the step; the vortex 
positions are disordered at the transitions and more ordered 
on the step. 
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FIG. 13: Vortex positions (filled dots), pinning site locations 
(open circles), and vortex trajectories (black lines) for the 
system in Fig. 12. (a) {tstep — t)/T — 0.06, below the locking 
step. The vortices are initially not locked to a specific angle, 
(b) In the center of the step region, {t — tstep) /T = 0.014, the 
vortices channel along 9 = 7r/3. (c) Above the locking step, 
{t — tstep) /T = 0.06, the vortex trajectories are disordered. 



step regions, the vortex lattice is almost completely tri- 
angularly ordered and contains only a small number of 
Zi — 5 and Zi = 7 paired defects. In Fig. 14(b) at 
t = tstep, the vortex lattice becomes strongly disordered 
and there is a proliferation of Zi — 5 and Zi = 7 defects. 
In the center of the velocity step at {t — tstep) /T = 0.014, 
Fig. 14(c) shows that the density of defects drops com- 
pared to Fig. 14(b). The remaining Zi — 5 and Zi — 7 
defects are aligned in the direction of the channeling 
motion, which can be seen from the trajectory plots in 
Fig. 13(b). The defect alignment indicates that on the 
steps, the vortex lattice has the characteristics of a mov- 



ing smectic. Above the step, as shown in Fig. 14(d) at 
(t — tstep) /T — 0.06, the vortices reorder into a nearly tri- 
angular lattice. A similar set of vortex lattice structures 
appear for each step and the topological defects align in 
the channeling direction on the steps. At the transitions 
on and off of each step, the vortex lattice disorders when 
a portion of the vortices move in a channeling direction 
while another portion of the vortices move in the driving 
direction. This combination of vortex motion disorders 
the lattice. 

We note that in general we do not observe complete 
locking, and even on a step, not all the vortices align 
with the drive as shown in Fig. 13(b). Much stronger 
velocity steps with complete locking appear for vortex 
motion in square pinning arrays at lower vortex densities 
of B/Bcf, — 1.05, as shown in Ref. [l^. For the system 
in Fig. 13, B/B^ = 2 and the locking effect is reduced 
since there are twice as many vortices as pinning sites. 
The vortex-vortex interactions tend to cause the vortex 
lattice to be triangular. The higher order velocity steps 
are small. In some cases, a high-order step disappears 
when the vortex- vortex interactions dominate to create a 
floating moving solid where the vortices form a triangular 
lattice and no locking effect occurs. In general we expect 
the steps to be more pronounced at lower vortex densi- 
ties. For higher vortex densities, only the velocity steps 
corresponding to the dominant symmetry directions of 
the pinning lattice will be discernible. 



B. Vortex Velocity Evolution with ac Amplitude 

As the ac amplitude increases, the velocity of the vor- 
tices increases and the locking step intervals diminish in 
width. Additionally, the number of steps that generate 
vortex disorder at their edges decreases, and for higher 
drives the disordering transitions occur only along the 
six prominent locking directions, as shown in Fig. 15 for 
A = 5.0. In Fig. 16 we compare the vortex lattice order- 
ing at A = 0.5 and A — 5.0 by simultaneously plotting 
Pg versus time for these two ac drive amplitudes. The 
change in lattice order when the vortices move along a 
symmetry direction of the pinning lattice is clear. At 
A — 0.5, the vortices are more ordered when moving 
along a symmetry direction, while at higher ac ampli- 
tudes this reverses and the vortices are less ordered when 
they move along a pinning symmetry direction. 

The overall change in Vx versus time for varied A is 
illustrated in Fig. 17. In Fig. 17(a) we plot Vx during 
the first half of the period for 0.16 < A < 0.43. One 
of the most striking features is that near t/T = 0.25 
for A = 0.35, which is highlighted by the darker curve, 
a new maximum in Vx develops that is associated with 
the random flow. This maximum becomes centered at 
t/T = 0.25 and increases in amplitude as A increases. 
An increasing fraction of the vortices depin during the 
random flow phase as A increases, and thus the magni- 
tude of Vx increases. The onset of the random phase is 



10 



(a) 



(c) 



(b) 



(d) 



FIG. 14: Voronoi constructions of vortex positions for the 
system in Fig. 13. White polygons are sixfold coordinated, 
Zi = 6; large black dots indicate polygons with = 5 and 
large gray dots denote polygons with Zi — 7. (a) Vortex con- 
figuration below the locking step at (tstep—t) /T — 0.06 at the 
point shown in Fig. 13(a). The vortex lattice is mostly tri- 
angular, (b) Vortex configuration at t = tstep, the transition 
into the locking phase. There is a proliferation of topologi- 
cal defects, (c) Vortex configuration at {t — tstep) /T ~ 0.014 
in the locking phase illustrated in Fig. 13(b). There is some 
partial ordering of the lattice, and many of the Zi = 5 and 
Zi = 7 defect pairs are aligned in the direction of the locked 
motion, (d) Vortex configuration at {t — tstep) /T = 0.06 above 
the locking step as shown in Fig. 13(c). The vortex lattice is 
ordered. 




FIG. 15: (a) Vx (dark line) and Vy (light line) versus time 
for the same system in Fig. 10 with A — 5.0. Here the step 
structures are very small, (b) The corresponding Pq versus 
time. Disordering transitions occur at the steps which are 
associated with the six symmetry angles of the pinning lattice, 
where some locking motion occurs. 



loidal assemblies when a colloidal crystal forms, it should 
be possible to observe the order-disorder phenomena we 
describe here. 



Commensurability Effects Near the B/B^ — 2.0 
State 



correlated with an increase in the level of fluctuations 
in Vx- For A > 0.435, there are no longer any locked 
sliding regimes; however, a small amount of locking per- 
sists and produces shoulders in Vx near t/T = 0.15 and 
t/T = 0.35 that gradually wash out at higher values of A. 
In Fig. 17(b) we plot K for 0.43 <A< 2.5, highlighting 
the transition from the random phase for the two lowest 
curves at ^ = 0.43 and A = 0.5 to the sliding phase at 
higher A where all the vortices are moving and new lock- 
ing phases for flow along the pinning sites occurs. This 
can be seen in the appearance and evolution of the lock- 
ing steps and the reduced fluctuations in Vx at higher 
A. 

Symmetry locking of the type shown in Fig. 10(b) for 
high drives has been observed previously both for vortices 
moving over periodic substrates [l^ and in experiments 
on colloids moving over a square substrate [5y|. In the 
colloidal studies, the dynamical locking occurred in the 
single particle limit, so the order-disorder type behavior 
we observe did not occur. We expect that for denser col- 



We next examine commensurability effects near the or- 
dered dimer state. In Fig. 18(a) we plot Vx for B/B^ — 
1.89 for a system with the same parameters as Fig. 1 
where we observed polarization switching of the dimer 
state. At the incommensurate filling, the large interstitial 
pins contain a mixture of dimers and monomers and the 
global orientational ordering of the dimers is lost. Those 
dimers that are present form some local ferromagnetic 
order in one of the three symmetry directions. Under 
an applied ac drive, a portion of the dimers align with 
the drive; however, the overall disorder in the sample cre- 
ates randomness in the energy barriers that a dimer must 
overcome to switch to a new orientation. This random- 
ness destroys the sharp global switching of the dimers 
shown in Fig. 3 at B/B^, — 2.0 and replaces it with a se- 
ries of smaller avalanche-like switching events. This is il- 
lustrated in Fig. 18(a) where only a portion of the dimers 
change orientation in any switching event. In Fig. 19(a) 
we show a time trace of the vortex trajectories for the 
system in Fig. 18(a) over several ac drive periods, illus- 
trating that the interstitial dimers and monomers remain 
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FIG. 16: Ps versus time for the system in Fig. 10 comparing 
the random phase at ^4 = 0.5 (lower curve), where the vortices 
are more ordered along the six symmetry directions, to the 
phase at A = 5.0 (upper curve), where the system is most 
disordered when sliding along the six symmetry directions of 
the pinning lattice. 



pinned but can move within each large interstitial region. 
The resulting trajectories are much more random than 
in the B/B^ = 2.0 switching case shown in Fig. 2. The 
switching pulses in Fig. 18(a) do not follow the same pat- 
tern for successive periods, indicating that the avalanche 
pattern is different from one period to another and has 
chaotic-like features. This type of chaotic avalanche be- 
havior is very similar to the Barkhausen noise effects 
found in disordered spin systems (st} . A fuller examina- 
tion of the Barkhausen or crackling noise effect observed 
in Fig. 18(a) will be presented in a later work. 

For B/B^ = 2.056, there is a mixture of dimers and 
trimers in the large interstitial regions, and for A = 0.075, 
a partial depinning of the interstitial vortices occurs cre- 
ating the large velocity bursts shown in Fig. 18(b). The 
bursts are generally confined to regimes where the parti- 
cles can channel along the six easy-flow directions of the 
pinning lattice, as illustrated in Fig. 19(b). For lower 
values of A, the partial depinning is lost and the system 
has a random switching behavior similar to that shown 
in Fig. 18(a). For very low ac amplitude, even the ran- 
dom switching is lost and the dimers and trimers remain 
locked in a fixed orientation. 



D. ac Effects on the Vortex Trimer State at 

B/B4, = 2.5 

As a function of field, higher order vortex molecular 
crystal states form whenever orientationally ordered ri- 
mer states occur, such as at B/B^ — 2.5, 3.0, or 3.5. 
These states can show a polarization effect as a func- 
tion of the external drive similar to that observed for 




FIG. 17: Vx versus time for different values of the ac ampli- 
tude showing the transitions from the random phases to the 
high ac drive phases, (a) From bottom to top, A = 0.16, 
0.225, 0.3, 0.325, 0.33, 0.35, 0.38, 0.405, and 0.43. The 
A = 0.35 curve is highlighed with a darker line, (b) From 
bottom to top, A = 0.43, 0.5, 0.6, 0.75, 1.0, 1.5, and 2.5. 



: 2.0. We focus on the case of B/B^ = 2.5, where 
the large interstitial regions each capture three vortices 
per pinning site and the trimers have a ferromagnetic or- 
dering pl^ . When the external drive is small, the trimers 
remain locked in one orientation without switching. For 
larger ac drives, switching of the trimer orientation oc- 
curs that is similar to the switching observed for the 
dimers. For even larger ac drives, a new type of switching 
effect occurs for the trimers that is absent for the dimer 
system. Within the large interstitial site, the three vor- 
tices in the trimer form a triangle. When the trimer is 
aligned along 9 = 7r/2 as in Fig. 5(c) of Ref. [i^ and 
the ac drive is also aligned along 6 = tt/2, the vortex at 
the top of the triangle experiences an extra force from 
the lower two vortices in the triangle and depins into the 
interstitial site above. The same process occurs at every 
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FIG. 18: (a) Vx vs time for a system with the same param- 
eters as in Fig. 2 with A = 0.075, but with B/B^ = 1.89. 
The vortices all remain pinned and the switching from one 
dimer orientation to another is random and reminiscent of 
Barkhausen noise, (b) versus time for the same system 
with B/B^ = 2.056, where a portion of the vortices begin to 
slide as the drive is rotated. 
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FIG. 20: Vy vs time for a system with the same parameters 
as in Fig. 2 but with B/B^ = 2.5, where a trimer state forms. 
From lower left to upper left, A = 0.075, 0.1, 0.125, 0.15, 
0.25, 0.3, and 0.4. For A < 0.15 there is a minimum in \Vy\ 
at t/T — 0.5, while for A > 0.25 a maximum in \Vy\ begins 
to form at t/T = 0. 
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FIG. 19: Vortex positions (filled dots), pinning site loca- 
tions (open circles), and vortex trajectories (lines) during 
several periods of motion, (a) The system in Fig. 18(a) at 
B/B^ = 1.89. The interstitial vortices follow disordered 
closed orbits, (b) The system in Fig. 18(b) at B/B^ = 2.056. 
The interstitial vortices depin and move from one intersti- 
tial region to another, while the vortices at the pinning sites 
remain pinned. 



large interstitial site and so there is an exchange of one 
vortex from each site to the site above. The net result 
is that the trimers reverse their orientation and are now 
pointing along 6 = 37r/2. The new orientation is more 
stable against depinning by the external drive. The two 
vortices at the top of the reoriented trimer produce a 
barrier for the motion of the vortex at the bottom of the 
trimer. At the same time, the top two vortices are too 
wide to hop out of the interstitial site into the interstitial 
site above. This can be regarded as a jamming effect. As 
the ac drive cycle continues, the drive rotates into the 
9 = 37r/2 direction, and the lower vortex in each trimer 



can jump into the interstitial site below, reorienting the 
trimer to point along 6 — 7r/2. Thus in this phase, trimer 
reorientation occurs by the exchange of individual vor- 
tices between neighboring interstitial sites, rather than 
by the rotation of the trimer itself. For higher ac drives, 
a sliding state appears when the vortices move along the 
symmetry directions of the pinning lattice. 

In Fig. 20 we plot Vy versus time for different values of 
A in a system at B / — 2.5 with the same parameters 
as in Fig. 1. At this filling a set of phases appear which 
are similar to those found at B/B^ = 2.0. At very low 
A, the trimers are pinned, while for larger A the system 
has sliding states interlaced with pinned states. A strong 
channeling effect produces local maxima and minima in 
Vy, similar to the B/B^ = 2.0 system. For A < 0.15, the 
vortices at the pinning sites remain pinned and a local 
minimum in \Vy\ occurs at t/T = 0.5. For A > 0.15 this 
minimum turns into a maximum when the vortices in 
the pinning sites begin to depin and generate a random 
flow phase. For A = 0.3, the transition from a locked 
flowing phase of interstitial vortices to the random phase 
produces a large jump in V^. In general, we find that 
as B/B^ increases, a more complicated set of features in 
the velocity force curves appears. 



IV. FIXED DC DRIVE IN THE y-DIRECTION 
WITH AN AC DRIVE IN THE x-DIRECTION 

In previous work on a honeycomb pinning array at 
B/B^ — 2.0, we found that when a dc drive is applied 
in the cc-direction, the vortices flow along either 9 — tt/& 
OT 9 = —it/6 depending on the original orientation of 
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the dimers in the ground state at F^; = [i^, \^ . Once 
the dimers are flowing under the dc drive, apphcation 
of an additional ac drive in the y direction can induce a 
switching between flow along 9 = Tr/6 and 9 ~ —tt/6 '4^. 
There is a threshold amplitude of the ac drive required to 
cause the switching in the symmetry breaking flow states. 
Here we consider the case of a dc drive F'^'^ = F'^'^y fixed 
in the y direction with an ac drive F"'^ = A sin(27rt/T)x. 
applied in the x direction. For A — 0, the dimers flow 
in one-dimensional paths along the y direction and there 
is no symmetry breaking of the flow, unlike the symme- 
try broken flow that occurs for an x direction dc drive. 
When an additional ac drive is applied in the x direction, 
transverse to the dc drive, the average value of Vy can be 
strongly influenced. 

In Fig. 21 we plot and Vy versus time for a system 
at B/B^ = 2.0 with a y direction dc drive of F''-'^ — 
0.15 and an x direction ac drive with A = 0.35. For 
< t/T < 0.075, the dimer motion is locked along 9 — 
tt/2. Although the magnitude of the ac drive is increasing 
in the positive x direction during this time period and 
therefore the overall net driving force vector is increasing 
in magnitude, Vy decreases monotonically with time until 
a switching event occurs at t/T — 0.075 when the dimers 
transition from flowing along = 7r/2 to flowing along 
9 = 7r/6. This switch coincides with a sharp jump in 
14 to a positive value and a cusp feature in Vy. After 
the switching event, Vx and Vy gradually increase with 
time until t/T = 0.25, when they begin to decrease. For 
0.075 <t/T < 0.4, the value of Vy is higher than it would 
be under only a dc force and no ac force, indicating that 
the application of a perpendicular ac drive can increase 
the velocity in the direction of the dc drive. For 0.25 < 
t/T < 0.5, both Vx and Vy decrease until Vx = Vy = 0, 
indicating that the system is pinned. The ac drive begins 
to increase in the negative x direction above t/T = 0.5, 
and as the dimers reorient to follow the net drive vector 
they realign with 9 = tt/2 and are once again able to slide 
along 9 = tt/2 for 0.525 < t/T < 0.58. At t/T = 0.58, 
another switching event occurs and the dimers slide along 
9 = ll7r/6. 

We term the pinning effect in Fig. 21 near t/T — 0.5 a 
jamming transistor effect. At t/T — 0.5, the ac force is 
zero in the a;-direction, so Vy would be expected to take 
the value it had at t/T = when the ac force was also 
zero. Instead, Fig. 21 shows that Vy = at t/T = 0.5. At 
this point, the ac drive has reoriented the dimers so that 
they are aligned along 9 — instead of along 9 = tt/2 
as they were at t/T = 0. In previous work on honey- 
comb pinning arrays at B/B^ = 2.0, we showed that if 
a dc drive is applied in the y-direction and the dimers 
are aligned perpendicular to the applied driving force, 
an unusually high critical depinning force appears due to 
the fact that the dimers are effectively jammed [13] ■ The 
pinned phase near t/T = 0.5 in Fig. 21 is an example of 
this jamming effect. The dimers act like rigid objects and 
cannot move through the y-direction channel between ad- 
jacent interstitial pinning sites unless they are aligned 
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FIG. 21: Vx (dark line) and Vy (light line) versus time for a 
system at B/B^ = 2.0 with a y direction dc drive of F'*'^ = 
0.15 and an x direction ac drive with A — 0.35. At t/T — 0, 
the dimers are moving along 6 — tt/2. At t/T = 0.075, a 
switching transition occurs and the dimers slide along 9 = 
tt/6. Near t/T = 0.5, a pinned phase (P) appears. This 
shows that the application of an ac drive induces a pinned 
phase in a system with a dc drive. We note that the ac drive 
amplitude is zero at t/T — 0.5 and the pinned phase occurs 
due to the reorientation of the dimers by the ac drive into a 
jammed configuration that is perpendicular to the dc drive. 
We call this the jamming transistor efi^ect. 



along 9 = tt/2. A key difference between the effect de- 
scribed here and the jamming in our previous work is 
that the jamming in Ref. occurred as an equilibrium 
transition in the ground state when the pinning force was 
increased. In the present case, the dimers are reoriented 
dynamically. In Fig. 22 we schematically illustrate the 
dynamic jamming effect. In Fig. 22(a), a dimer oriented 
along 9 = tt/2 can move between the pinning sites and 
has a low depinning threshold. If the dimer is oriented 
along 6* = as in Fig. 22(b), parallel to the two pinning 
sites, the energy barrier for the dimer to pass between 
the pinning sites is much higher and the depinning force 
is correspondingly higher. 

In Fig. 23 we plot Vx and Vy for a system contain- 
ing only interstitial monomers at B/B^ = 1.5 under the 
same drive configuration as in Fig. 21, with a dc drive 
applied in the y-direction and an ac drive applied in the 
x-direction. There is a modulation in Vy, with a dip near 
t/T = 0.12 at the point where Vx rises above zero; how- 
ever, Vy never drops to zero and there is no pinned phase. 
At t/T — 0.5, Vy reaches the same value it had at t/T — 
at the beginning of the period when the ac drive compo- 
nent was zero. This is distinct from the behavior of the 
dimer system, where Vy dropped to zero at t/T — 0.5 due 
to the reorientation of the dimers. We find behavior sim- 
ilar to that shown in Fig. 23 over a wide range of ac and 
dc drive parameter values for B/B^ = 1.5. In general, 
the application of a force in the x-direction to a system of 
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(a) (b) 

FIG. 22: Schematic of the jamming effect. Light circles: 
dimer vortices in interstitial site. Dark circles: pinned vor- 
tices. Open circles: pinning sites, (a) A dimer oriented along 
6 = TT /2 can easily move between two occupied pinning sites 
and has a low depinning threshold, (b) A dimer oriented 
along 9 = has difficulty moving between the occupied pin- 
ning sites, and the depinning threshold is much higher. 
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FIG. 23: Vx (dark line) and Vy (light line) for the same system 
in Fig. 21 with B/B^ = 1.5 so that no dimers are present. The 
dc drive is in the y direction with F'''^ = 0.425 and a sinusoidal 
ac drive is applied in the x direction with A = 0.35. Unlike 
the case of B/B^ = 2.0, the velocity responses at t/T = 
and t/T = 0.5 are the same. 



monomers being driven in the y direction does not induce 
a pinned phase or a phase where Vy = Q. The fluctua- 
tions in Vy increase when Vx becomes finite because the 
flow in this regime is disordered, although the vortices 
at the pinning sites remain pinned. The monomers ran- 
domly alternate which side of the occupied pinning sites 
they move around. 



A. Jamming Transistor Effect 

A system in which a dc velocity component can be 
switched on and off with an ac drive applied perpendic- 
ular to the dc drive has similarities to field effect tran- 
sistor devices. The strong cusp in Vy near t/T — 0.075 
in Fig. 21 and the corresponding sudden jump in Vx sug- 
gest that the rapid switching time could be used to create 
semiconductor type devices similar to the transverse field 



effect depinning threshold modulation in charge density 
wave systems O E^- Since the external ac field can 
control the dimer orientation and since the depinning 
threshold depends on the orientation of the dimers, the 
system should exhibit strong memory effects of its origi- 
nal configuration. For example, different settings of the 
initial ac and dc drive values would orient the dimers in 
a particular direction and the pinning threshold would 
depend on the history of the preparation of the dimer 
states. These field effect type properties could be use- 
ful for making vortex-based devices or could be used as 
models for constructing similar systems on smaller scales 
using ions or Wigner crystals. 

We showed in Fig. 23 that monomers do not exhibit 
a jamming effect. In general, the jamming effect occurs 
when B/Bfj, is close to B/B^ = 2.0. As long as there 
are enough interstitial dimers to percolate throughout 
the entire sample, the dimers can induce a jamming ef- 
fect which has a side effect of pinning any monomers or 
trimers that are present. Even at B/B^ — 2.0, there is 
a threshold value of F'^'^ above which the jamming tran- 
sition no longer occurs, as we describe in Section IV C. 
For fields away from B/B^ = 2.0, this dc threshold force 
is reduced. For B/B^ < 1.75, a complete jamming effect 
is lost; however, portions of the system populated with 
dimers can still jam. For B/B^ < 2.0 and low enough 
F'^'^, the monomers are simply pinned at the interstitial 
sites and jamming of the dimers has no effect on the 
monomer pinning. 



B. Varied ac Amplitude and Fixed dc Amplitude 

We next examine the effect of varied ac amplitude for a 
system with fixed F'^'^ = 0.1 applied in the y direction at 
B/B^ = 2.0. If the ac amplitude is too small, the vortex 
motion remains locked along = 7r/2 in the y-direction. 
In Fig. 24(a) we plot Vx and Vy for A = 0.16. At this 
ac drive value, the jump in Vx is close to its maximum 
amplitude, and the corresponding cusp feature in Vy is 
reduced in size. The pinned phase near t/T = 0.5 is also 
more extended than in the sample illustrated in Fig. 21. 
In Fig. 24(b) we show 14 and Vy in the same system at 
A = 0.15, which is just below the ac threshold required 
to unlock the motion from the y-direction, = 11/2. In 
this case Vx = throughout the cycle while Vy has a 
periodic modulation but is always finite. Even though 
the vortices are not translating in the x-direction, the 
application of an ac force in the x-direction causes the 
vortices flowing along the y direction to shift toward one 
side of the 6 — it/ 2 flow channel. This increases the 
strength of the interaction between the moving vortices 
and the pinned vortices and increases the effective drag 
on the moving vortices, reducing Vy whenever the abso- 
lute value of the ac drive is maximum in the positive or 
negative x-direction. In Fig. 25 we show a series of Vy 
versus time curves for increasing ac drive amplitude with 
A = 0.35, 0.45, 0.55, and 0.85. As A increases, the tran- 
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FIG. 24: (a) 14 (dark line) and Vy (light line) vs time for 
a system with a y direction dc force of F'^'^ — 0.15 and an x 
direction ac force with A — 0.16. Here the jamming transition 
appears along with a strong switching effect in Vx- (b) The 
same system at A = 0.15, where the ac amplitude is too small 
to induce a switching effect and Vx = 0. 
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FIG. 25: Vy vs time for the same system in Fig. 21 with 
F'^'' = 0.15 at ^ = 0.35, 0.45, 0.55, and 0.85, from bottom to 
top. The large jump for A > 0.35 near t/T = 0.25 is the onset 
of the random phase where vortices at the pinning states start 
to depin. 



grows until for A < 0.15 the system is always locked 
to 6* tt/2, as shown in Fig. 24(b) for A = 0.15. For 
0.15 < A < 0.45, the system passes through 9 — Tr/2 
and 9 — IT / 6 locked flow phases and a jammed or pinned 
phase in half a period. The random flow phase flrst ap- 
pears for A > 0.45. In the random flow phase, vortices at 
the pinning sites depin when the net driving force from 
the combined dc and ac drives is maximum, which occurs 
at t/T = 0.25. The jammed phase which is centered at 
t/T = 0.5 becomes narrower but still persists for increas- 
ing A. For A > 1.25 (not shown) there are new types 
of moving phases that arise within the random phase re- 
gion when all the vortices depin and start sliding along 
the pinning sites, producing the order-disorder transi- 
tions studied earlier. 



sition from = 7r/2to0 = 7r/6 motion occurs at earlier 
t/T. For A > 0.35, there is a large jump into a fluctu- 
ating phase centered near t/T — 0.25, corresponding to 
the depinning of a portion of the pinned vortices into a 
random flow regime. The system can then reorganize to 
the 9 = Tr/6 flow phase at a later point in the period 
when the force from the ac drive goes to zero and the net 
driving force is no longer large enough to depin vortices 
from the pinning sites. The vortices enter a pinned phase 
at t/T = 0.5. 

In Fig. 26 we show the dynamical phase diagram of 
time versus ac amplitude A during half of the ac drive 
period with flxed F'^" = 0.15. At small t/T when the 
ac force is first increasing from zero, the system always 
starts with the dimers flowing along 9 = 7r/2. For de- 
creasing A, the temporal extent of the 9 = tt/2 region 



C. Constant ac Amplitude and Varied dc 
Amplitude 

We next consider the case where the x direction ac 
drive amplitude is flxed at A = 0.2 and the y direction dc 
drive amplitude is varied at B/B^ = 2.0. At F'^'^ = 0, the 
behavior of the system resembles that observed in earlier 
studies with only a dc drive in the x-direction [4^, [4^ . 
In this limit the ac drive applied along the x direction in- 
duces flow along 9 — tt/Q or 9 — llvr/G, with the flow re- 
versing during the second half of the ac cycle to 9 = 77r/6 
or 9 — 5tt/6, respectively. In Fig. 27(a) we plot Vx and 
Vy for a system with F'^'^ = and A = 0.2. Here the vor- 
tices are initially pinned and then undergo a transition 
to the moving 9 = lln/Q phase. In this case, the dimers 
were oriented along 9 = ll7r/6 in the initial pinned 
ground state. The amplitude of Vx is larger than Vy since 
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FIG. 26: The dynamical phase diagram for time tjT vs ac 
amplitude A for the system in Fig. 21 with a fixed dc drive of 
F^'^ — 0.15 in the y-direction and an ac drive applied in the 
x-direction. In the 9 = 11/2 phase the dimers flow along the 
y-direction, while for higher A they flow along 6 — n/6. The 
pinned phase is marked P while R is the random phase where 
the vortices start to depin from the pinning sites. 



|14| = |Acos(ll7r/6)| and \Vy\ = |Asin(ll7r/6)|. When 
the ac drive reverses, the dimers flow along 9 — 5tt/Q. 
The symmetry breaking of the dimers in the ground state 
fixes the flow direction of the dimers for all later times. 
If we consider an incommensurate case where there are 
additional disordered phases, the system dynamically or- 
ganizes into one of the broken symmetry phases. In this 
case, when the ac drive reverses sign, the system passes 
though the disordered phase which destroys the mem- 
ory of the dynamically broken symmetry phase, and the 
vortices have an equal chance of dynamically breaking 
symmetry in either of the two possible directions. 

At B/B^ = 2.0, if a small but finite dc drive F'^" is 
applied in the y direction, it is unable to break the sym- 
metry of the ground state, and the ground state dimer 
orientation determines the resulting dynamics, as illus- 
trated in Fig. 27(b) for a system with F'^'^ = 0.05. Here 
the vortices are initially pinned and then depin to fiow 
along 9 — 7r/6, reflecting the orientation of the dimers 
in the pinned state. During the second half of the drive 
period, the dimers fiow along Ttt/G and \Vx\ and \ Vy \ are 
both lower than they were during the first half of the 
drive period. This is due to the biasing force of the dc 
drive component which is applied in the positive y di- 
rection. When F'^'^ is sufficiently strong, the negative 
component of the y velocity during the second half of the 
drive period is lost. 

In Fig. 28 we plot Vy versus time for A = 0.2 and F'^'^ — 
0.05, 0.125, 0.2, 0.225, and 0.275. For F'^^ > 0.125, Vy 
is always positive. For F'^'^ > 0.1, the vortices initially 
move along 9 = tt/2 and as F'^'^ increases, the magnitude 
of Vy increases, as does the magnitude of the jump from 
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FIG. 27: Vx (dark line) and Vy (light line) vs time for a 
system at B/B^ = 2.0 with a y direction dc drive and an x 
direction sinusoidal ac drive with A — 0.2. (a) F'^'^ — 0. Here 
the vortex motion is locked along S = ll7r/6 or 57r/6 during 
a portion of the first and second halves of each ac drive cycle, 
(b) F'^'^ — 0.05. Here, the dc drive is not strong enough 
to lock the motion to only the positive y-direction, and the 
direction of motion is determined by the symmetry breaking 
of the dimers in the ground state. 



9 = 7r/6 back to 6* = 7r/2 near t/T = 0.5. At F'^" = 
0.275, the 9 = 7r/6 phase is replaced with a randomly 
fluctuating phase when the net force from the combined 
dc and ac drives is strong enough to depin vortices from 
the pinning sites, similar to the behavior shown earlier. 
For F'^'' < 0.2, Vx (not shown) steadily rises with t/T 
and reaches a maximum value near t/T = 0.25 which 
saturates for F'^'^ > 0.2. The saturation occurs since the 
ac component in the x-direction is flxed. The onset of the 
random phase also produces increased fluctuations in Vx 
at F'^" = 0.275. 

In Fig. 29(a) we show the continued evolution of Vy 
from Fig. 28(a) for F"^" = 0.33, 0.41, 0.45, 0.5, 0.55, and 
0.6. The = 7r/2 sliding phase diminishes in size with 
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FIG. 28: Vy vs time for the same system in Fig. 27 with an 
X direction ac drive of ^ = 0.2 and a y direction dc force. 
From bottom to top, F'^" = 0.05, 0.125, 0.2, 0.225, and 0.275. 
The onset of the random phase where vortices depin from the 
pinning sites appears as a jump into a strongly fluctuating 
phase for the F"^" = 0.275 curve. 



increasing F'*'^ and vanishes completely for F''^'^ > 0.475 
while the random flow phases grow in extent. Within 
the random flow phases for F^c > 0.475, there is a mod- 
ulation of Vy with dips near t/T = 0, t/T = 0.5, and 
t/T = 1, along with smaller modulations near t/T — 0.25 
and t/T = 0.75. It is interesting to note that the ratio 
of Vy at t/T = 0.25 to Vy at t/T = 0.5 reaches a maxi- 
mum at F''-^ = 0.45. Here Vy is strongly reduced when 
the system enters the 9 = n/2 phase because the vor- 
tices that depinned during the random flow phase repin 
in the 9 = tt/2 phase and no longer contribute to V^^. In 
Fig. 29(b) we show the corresponding Vx curves which 
indicate that the 9 = tt/2 locked flow phase disappears 
for F'^" > 0.475. 

In Fig. 30 we plot the dynamic phase diagram for 
F'^'^ versus time at fixed A = 0.2 for the system in 
Figs. 28 and Fig. 29 at B/B^ = 2.0. For t/T < 0.1 
and F'^'^ < 0.075, the vortices are initially in the pinned 
phase. As t/T increases and the ac force becomes larger, 
the vortices organize into a pinned symmetry broken 
state for F'''^ < 0.075 or dcpin into the 9 = tt/Q sliding 
state for 0.075 < F'^" < 0.475. The initial pinned state is 
lost for F'^'^ > 0.1 when the dc drive is strong enough to 
cause the dimers to slide along 9 = ■7r/2. For F''-'^ > 0.475, 
the dc drive combined with the motion of the interstitial 
vortices is strong enough to dcpin the vortices at the pin- 
ning sites and the system enters the random phase even 
when F""" = at t/T = 0. For 0.1 < F^"" < 0.25, the 
system switches from = 7r/2 to = 7r/6 dimcr flow, 
while the 9 = tt/6 phase is lost for F'^'^ > 0.25 and the 
system enters the random flow state. The second pinned 
state near t/T = 0.5 disappears for F'''^ > 0.2. Due to 
the polarization of the dimers, this is higher than the 
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FIG. 29: (a) Vy vs time for the system in Fig. 28 at ^ = 0.2 
and = 0.33, 0.41, 0.45, 0.5, 0.55, and 0.6 (from bottom to 
top). For F'^^ > 0.45, the shding phases 9 = ti/2 at t/T = 
and t/T = 0.5 are lost, (b) The corresponding 14 vs time at 
F"^" = 0.33, 0.41, 0.45, 0.5, 0.55, and 0.6 (from bottom to top 
at t/T = 0). 



drive at which the pinned state disappears for t/T = 0. 
For higher dc amplitudes than those shown in Fig. 30, all 
of the vortices depin and slide along certain directions of 
the pinning lattice, forming the order-disorder transitions 
studied earlier. 



D. Commensurability Effects 

We next consider the effects of varied B/B^. Figure 31 
shows Vy versus time for B/Bs = 2.44, 3.0, 3.4, and 
4.72 for a system with fixed F"''^ = 0.35 applied in the 
y-direction and an ac force with A = 0.2 applied in the x- 
direction. For clarity, the velocities are normalized by the 
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FIG. 30: Dynamic phase diagram F'*'^ versus time for the 
same system in Figs. 28 and Fig. 29 at B/B^ = 2.0 under 
an X direction ac drive with A — 0.2 and a y direction dc 
drive. P: pinned phase; SB: symmetry breaking phase where 
the symmetry is broken in the pinned dimer configuration; 
R: random phase where vortices at the pinning sites become 
depinned. There are two sliding phases along 9 = 7r/2 and 
6 = Tv/6. 
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FIG. 31: Vy vs time at varied B/B^ with a y direction dc 
drive of F"^" = 0.35 and an X direction ac drive with A = 0.2. 
For clarity, the velocities are normalized by the value of Vy 
at B/B^, = 2.0. From top to bottom, B/B^, = 2.44, 3.0, 3.4, 
and 4.72. For B/B^ — 3.4, a complex cusp structure forms 
and two different types of ordered phases appear. 



value of Vy at B/B^ = 2. For B/B^ = 2.44, the system 
is in the random phase throughout the driving period; 
however, a prominent channehng effect persists as shown 
by the minima at t/T = 0.5 and t/T = 1.0 as well as 
the smaller shoulders at t/T = 0.25 and t/T = 0.75. For 
B/Bij, = 3.0, the system forms several new ordered flow 
phases visible as regions of reduced fluctuations in Vy. 
The ordered interstitial flow phase ends at t/T — 0.04 
with a transition into a random flow regime, as shown by 
the drop in Vy. The same ordered interstitial flow state 
reappears just before and after t/T — 0.5, as indicated by 
the sharp cusps in Vy, and also occurs for B / B^ = 3.4. In 
Fig. 32 we illustrate the vortex trajectories in the ordered 
interstitial flow phase for B/B^ — 3.4. The moving in- 
terstitial vortices channel in periodic winding paths at an 
angle between the occupied pinning sites. At this partic- 
ular fllling, a small portion of the interstitial vortices are 
trapped behind the occupied pinning sites. The ordered 
interstitial flow phase at B/B^ — 3.0 is nearly identi- 
cal except that the immobile interstitial vortices are not 
present. 

In Fig. 31, as the ac force increases for B/B^ = 3.0 the 
vortices at the pinning sites depin into a floating triangu- 
lar structure and the flow becomes random at t/T = 0.1. 
Near t/T — 0.15, the fluctuations in Vy are reduced when 
a portion of the vortices lock into a channeling motion 
along the pinning rows. For B/B^ = 3.4, the same float- 
ing triangular lattice appears; however, the transition 
from the random phase to the channeling phase centered 
at t/T — 0.25 is much sharper. The plot of Pg versus 
time for the B/B^ = 3.4 system in Fig. 33 shows the 
transition to the floating triangular state, and also indi- 




FIG. 32: Vortex positions (filled dots), pinning site locations 
(open circles), and vortex trajectories (lines) in the ordered 
interstitial flow phase at t/T = 0.025 for the B/B4, = 3.4 
system in Fig. 31. Here the moving interstitial vortices form 
ordered channels between the pinning sites. A small portion 
of the interstitial vortices are pinned behind occupied pinning 
sites. This same type of ordered interstitial fiow phase occurs 
at B/B^ = 3.0. 
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FIG. 33: The fraction of six- fold coordinated vortices Pq vs 
time for the B/B^ = 3.4 system from Fig. 31. The cusp 
structures in Vy in Fig. 31 correlate with the features in Pq. 
In the channeling phase centered at t/T = 0.25, the vortices 
form a mostly triangular lattice. 



cates that the vortex lattice regains some order during 
the ordered interstitial flow phases. The transitions into 
the disordered regimes appear as drops in Pg. 

For B/B^ = 4.72 in Fig. 31, random flow phases dom- 
inate the behavior; however, several cusp features still 
appear in Vy . The vortex lattice is also mostly triangular 
near t/T = 0.25. A new feature at this field is the trans- 
formation of the minima in Vy at t/T — 0.5 into a small 
maximum. These results show that there is a complex 
variety of phases at the different fields which have clear 
transport and ordering signatures. 



V. DISCUSSION 

In this work we considered varied ac drive amplitude 
but kept the ac drive frequency fixed. Our results should 
be robust in the low frequency regime, and for frequencies 
lower than the frequency considered here, the results do 
not change. At much higher frequencies, the system does 
not have time to respond to the ac drive and the switch- 
ing events will be lost. In our previous work, we found 
that for a frequency about 10 to 100 times higher than 
the frequency considered here, the system stops respond- 
ing to the ac drive [i^ . In Fig. 34 we show an example of 
this effect for a system with a dc drive applied in the y- 
direction and an ac drive applied in the a;-direction, with 
the same parameters as the system in Fig. 23 but at a 
higher frequency. For low frequency, the system shows a 
switching effect from 9 = tt/2 flow to 6 = tt/6 flow. In 
Fig. 34 the frequency is 10'^ times higher and the vortex 
trajectories show oscillations; however, the net motion 
remains locked to 9 = 7r/2. 

Our previous work with this system showed that tem- 
perature effects can change the phase boundaries slightly. 
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FIG. 34: Vortex positions (filled dots), pinning site locations 
(open circles), and vortex trajectories (lines) during several 
drive periods for the same system in Fig. 23 with a y direction 
dc drive of F"^" = 0.425 and an X direction ac drive with 
A = 0.35. In Fig. 23 the frequency was 10'^ times lower and 
the system exhibited a switching from 6 = Tr/2 flow to 9 = n/6 
flow. Here, the ac frequency is higher than the ability of 
the vortices to respond, and the motion stays locked in the 
9 — Tv/2 direction. 



but that the ordered dynamical phases persist until the 
melt ing transition of the interstitial vortices is reached 
[H, SSl- In the present work, we focused on the limit 
where the pinning is strong enough that there is a regime 
in which the vortices at the pinning sites can remain 
pinned while the interstitial vortices flow. As Fp de- 
creases, many of the dynamical phases will disappear as 
the dimers lose their independent effective orientational 
degree of freedom. For low enough Fp, additional transi- 
tions can occur in the vortex lattice where a portion of the 
formerly pinned vortices are no longer located at the pin- 
ning sites, producing a triangular structure [3^. [45l [ssj. 



VI. SUMMARY 

We have shown that in honeycomb pinning arrays, 
vortex molecular crystalline states which have effective 
orientational degrees of freedom can have their orienta- 
tion controlled by an external ac drive. We speciflcally 
focused on vortex dimers at the second matching fleld 
which have three different possible orientations. When 
an ac drive is applied in both the x and y directions, the 
orientation of the dimers can be made to follow the ex- 
ternal drive while the vortices remain pinned in the large 
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interstitial site. We term this a polarization effect. In this 
regime the dimcrs rotate by a combination of continuous 
motion and a switching type behavior that appears as 
large velocity jumps when the dimers switch to the dif- 
ferent symmetry angles of the pinning array. For very 
weak ac drive amplitudes, the dimers remain locked in 
their initial orientation, showing that there is a thresh- 
old force needed to induce the pinned switching behav- 
ior. For large ac drive amplitudes, the dimers depin and 
flow between the pinned vortices along certain symmetry 
directions of the pinning lattice, producing distinctive 
features in the x and y components of the velocity dur- 
ing each ac period. For even larger ac drive amplitudes, 
all the vortices depin and the system exhibits a series of 
order-disorder transitions and symmetry locking effects 
as a function of ac drive. We also consider the case of 
a fixed dc drive applied in one direction and an ac drive 
applied in the perpendicular direction. Here, the applica- 
tion of the ac drive can cause the moving dimers to orient 
in such a way that they become effectively jammed. We 



call this effect a jamming transistor and show that it oc- 
curs due to the fact that the dimers can slide more easily 
if they are oriented along the symmetry axis of the pin- 
ning rows. The ac drive can also induce a pronounced 
dynamical switching that can be seen in both velocity 
components. Many of the switching features may be use- 
ful for creating new types of vortex based devices or flux- 
tronics, and can serve as examples for other systems in 
which molecular crystal type states can form, such as 
ions or Wigner crystals on periodic substrates, colloidal 
molecular crystals, or vortices in Bose-Einstein conden- 
sates. 
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